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Abstract. We describe an algorithm to count the number of distinct real zeros of 
a polynomial (square) system /. The algorithm performs C'(log(nD«(/))) iterations 
(grid refinements) where n is the number of polynomials (as well as the dimension of 
the ambient space), D is a bound on the polynomials' degree, and k(/) is a condition 
number for the system. Each iteration uses an exponential number of operations. The 
algorithm uses finite-precision arithmetic and a major feature in our results is a bound 
for the precision required to ensure the returned output is correct which is polynomial 
in n and D and logarithmic in k(/). The algorithm parallelizes well in the sense that 
each iteration can be computed in parallel time polynomial in n, logD and log(«;(/)). 

1 Introduction 

In recent years considerable attention was put on the complexity of counting problems over 
the reals. The counting complexity class #Pu was introduced [501 and completeness results 
for #Pr were established [3 for natural geometric problems notably, for the computation 
of the Euler characteristic of semialgebraic sets. As one could expect, the "basic" ^Pjf- 
complete problem consists of counting the real zeros of a system of polynomial equations. 

Algorithms for counting real zeros have existed since long. One such algorithm follows 
from the work of Tarski [25| on quantifier elimination for the theory of the reals. Its complex- 
ity is hyperexponential. Algorithms with improved complexity (doubly exponential) were 
devised in the 70s by Collins |5j and Wiitrich [27]. A breakthrough was reached a decade 
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later with the introduction of the critical points method by Grigoriev and Vorobjov [T31[T2] 
which uses exponential time. Algorithms counting connected components (and hence, in 
the zero-dimensional case, solutions) based on this method can be found in [TH [TB] , and in 
the straight-line program model of computation in |T] . These algorithms parallelise well in 
the sense that one can devise versions of them working in parallel polynomial time when 
an exponential number of processors is available. The ^P^-completeness of the problem 
strongly indicates that this is the best we can hope for. 

All the algorithms mentioned above are "symbolic algorithms." They have been devised 
upon the premise that no perturbation or round-off error is present. Were this not the 
case, it is not difficult to see that errors would accumulate quite badly. Roughly speaking, 
these algorithms construct some object of exponential size on which some basic computation 
(e.g., linear algebra) is eventually performed. A question is posed, can one devise "numerical 
algorithms" (maybe iterative, which need not terminate for ill-posed inputs) with a better 
behavior viz the accumulation of round-off errors? For the problem of deciding the existence 
of (or computing) a zero of a polynomial system such algorithms were given in [HI [6l I18j . 
The goal of this article is to describe and analyze a numerical algorithm for zero counting. 
We will do so by developping appropriate versions of the tools used in [SI |2] . 

Let di, . . . ,dn £ N and d = (di, . . . , c?„). We will denote by Tid the space of polynomial 
systems / = (/i, . . . , /n) with fi e M[A'o, . . . , Xn] homogeneous of degree di. 

Zero rays of polynomial systems / G Hd are associated to pairs of zeros {—(, C) of the 
restriction f^gn of / to the n-dimensional unit sphere S"" C M""*"^. Thus, it will be convenient 
to consider a system / G Hd as a (central symmetric, analytic) mapping of S"" into M". If 
we denote by Z{f) = S 5" : f{Q = 0} the zero-set of / in 5" then the number #r(/) of 
zero rays of the system / is half the cardinality of Z{f). 

In this paper we describe a finite-precision algorithm computing #r(/), given / G Hd- 
To analyze its complexity and accuracy, besides the number n of polynomials, we will rely on 
two more additional parameters. One is D = maxi<„ di. The other is a condition measure 
k(/) for the system /. We will describe this measure in detail in Section [2] below. We 
will also let S = max Si where 5*^ is the number of non-zero coefficients of fi. Note that 
S is bounded by a simple expression in terms of n and D, namely, S = ("^'^). Yet, we 
will express dependancy on S since this may be relevant for the case of sparse systems of 
polynomials. Our main result is the following. 

Theorem 1.1. There exists an iterative algorithm which, with input f G Tid, 

(1) Returns #r(/). 

(2) Performs C'(log(nDK(/))) iterations and has a total cost (number of arithmetic opera- 

tions) of 

o(,„g(„D«(/))(„ + i,^(H(ii±ii5W)^"), 

where a* ~ 0.0384629388 . . . is a universal constant. 

(3) Can be well-parallelized in the sense that it admits a parallel version running in time 

0{n^ ln(nDK(/))(ln(nDK(/))2 + \n{a,f)) 
with a number of processors exponential in this quantity. 
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(4) Can be implemented with finite precision (both versions, sequential and parallel). The 

running time remains the same (with a* replaced by a, ~ 0.028268 ■ ■ ■ ) and the re- 
turned value is #r(/) as long as the machine precision (i.e., the round-off unit) u 
satisfies 

1 

u < 

- O (D2n5/2K(/)3(logS' + n3/2D2K(/)2)) ■ 

(5) It can be modified to return, in addition and for each real zero ( e S*" of f , an approx- 

imate zero X of f in the sense that Newton's iteration, starting at x, converges to C 
quadratically fast. 

Remark 1.2. (i) A system / for which arbitrarily small perturbations may change 
the value #k(/) is considered ill-posed in our context since for arbitrarily small machine 
precisions finite precision algorithms may return an incorrect value. Consequently, the 
condition number «;(/) is infinite in these cases (and only then). This happens when / has 
multiple real zeros and, in particular, when / has infinitely many real zeros. In these cases 
the algorithm of Theorem 1 1 . 1 1 mav not halt. 

(ii) Numerical algorithms compute functions ip on real data. Error analysis for algorithms 
computing (vectors of) real numbers — i.e. for which the image of f has non-empty interior — 
are usually expressed in terms of bounds for the relative error of the computed quantities. 
That is, for data d, bounds in 

Md)~fiiipid)\\ 
Md)\\ 

where fl{ip{d)) is the vector actually computed with finite precision. This relative error 
varies continuously with d and depends on the condition of d and on the precision u. Such a 
form of analysis, however, becomes meaningless when computing quantities taking a finite 
number of values. Indeed, if Ra denotes the set of input data d for which ip{d) = a the 
following happens. When d is in the interior of Ra we have that the relative error above 
is for sufficiently small u. In contrast, when d is on the boundary of Ra, that error may 
remain constant for all u > 0. Because of this, error analysis for this kind of discrete-valued 
problems has a different form, as in Theorem ll.il One bounds how small u needs to be to 
guarantee a correct answer. Such a bound, needless to say, also depends on the condition 
of the data d. Examples of this type of analysis can be found in [H [SI [71 [5] . In each of these 
references a condition number for the problem at hand occurs in the error analysis. We note 
that the one in [6] is essentially our k(/). 

The rest of the paper is organized as follows. In Section[2]we describe the basic objects we 
will deal with as well as fixing the notation. In Sections [3| and [4| we prove the two technical 
results our algorithm relies on. In Sectior(5|we describe the algorithm under the assumption 
of infinite precision and we prove parts (1), (2), and (3) of Theorem 11.11 The geometric 
ideas making the algorithm work are best seen in this context. Section [5| then describes the 
necessary modifications to make the algorithm work as well under finite precision. These 
modifications are simple and can be summarized by saying that we relax a bit the inequalities 
tested in the algorithms to make room for the finite-precision errors to fit in. 

2 Preliminaries 

Denote by Hd the subspace of ^[Xo, . . . , X„] of homogeneous polynomials of degree d. Then, 
Ti-d = Ti-di X ■ ■ ■ X "Hd^ ■ 
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If g G Ti-d we write 



where J = (Jo, . . . 



, Jn) is assumed to range over all multi-indices such that \J\ = X]fc=o ~ 
. . and gj G R. Multinomial coefficients are defined by: 



dl 



JolJll---Jnl' 

The space Tid is endowed with the inner product 

gjhj 

\J\= 



{9,h) 



rd\ 
\=d \J) 



which gives rise to the norm ||g|| = {g, g). These norms, for di, . 
Hd by taking for / = (/i, . . . , /„) e Hd- 



. ,dn, induce a norm in 



11/11 -Il(/i 



I fn 



max II/, 

l<2<n 



Let 0{n + 1) be the orthogonal group. The inner product above is known to be 0{n + 1)- 
invariant: for all Q e 0(n + 1) and all 5, /i G Hd, 

{goQ,hoQ) ^ {g,h). 

(This is a direct consequence of |26[ IIT7] or [3J Theorem 1 p. 218], by considering 0{7i + 1) 
as subgroup of U{n + 1)). The associated norm ||/|j on Hd is therefore also 0{n + 1)- 
invariant. We will use this norm on Hd all along this paper. For x = {xi, . . . , Xn) G K" we 
recall that ||a::||2 = (xf + • • • + x^y/^ and ||a;||oo = maxdxil, . . . , \xn\}- We will often denote 
||x||2 simply by ||a::||. 

For / G and X G S*" define 



.(/,a:) = 11/11 



(1) 



where Df{x)\T^S" is the restriction to the tangent space of x at S" of the derivative of / at 
X and the norm is the spectral norm, i.e. the operator norm with respect to || II2. We now 
define the condition number K{f) of / G Hd- 



max mm < ^norm 



l/ll 



1/(2^)11 



Remark 2.1. The quantity k(/) is closely related to other condition numbers for similar 
problems. 

A version of the quantity /inorm(/, C) was introduced in [2Tl[22j[23j (see also [H Chap- 
ter 12]) for a complex polynomial system / and a zero C of / in the complex unit sphere 
5^ C C"^^. The normalized condition number of such a system / was then defined to be 



.(/) 



niRX /^norm 

CeS^I/(C)=o 



(/,C). 



(2) 
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Actually, the version of /Xnorm(/, C) introduced in [5T1 [551 US] differs from ([T]) in the fact 
that 11/11 is defined as £ ll/if )^^^ (and there is no ^/n factor) . It is bounded above by the 
expression in ([T]). 

Over the reals, the right-hand side in ^ may not be well-defined since the zero set 
of / may be empty. In [8] real systems were considered (as in the present paper) and 
an algorithm deciding feasibility of / (i.e., whether / has a real zero) was proposed. Its 
complexity was analyzed in terms of a condition number which, using our notation and 
modulo minor details, is defined as follows 

min /inorm ( f , C) if f is feasible 
li/li 

max ,, — if f is infeasible. 

I Ces" |l/(C)||oo 

Note the use of min (instead of max) in the first line above. This is due to the fact that the 
time needed for the algorithm in [5] to detect the existence of a zero depends on the best 
conditioned zero of /. The existence of other, poorly conditioned (or even singular), zeros 
of / is irrelevant. 

Shortly after, the algorithm in [5] was extended to an algorithm which would, in addition 
and if / is feasible, return a zero of / [6]. The complexity of this extension was studied 
in terms of a condition number (denoted gi(/) in [S]) which, essentially, coincides with our 
«(/)■ 

Proposition 2.2. For all f e Hd, «(/) > 1- 

Proof. Let x E S"^. Because of orthogonal invariance, we may assume without loss of 
generality that x — cq :— (1, 0, . . . , 0). 

It is then immediate that |j/(x)||oo < ||/||- This shows that the second expression in the 
definition of k is at least 1. 

For the first expression, i.e., Ainorm(/, a:), define g = {gi, . . ..g^) g TLd hy gi{X) = fi(X)~ 
fi{eo)XQ\ Then ^(eo) — and [H Corollary 3 p. 234], /Xnorm(ff,eo) > 1 (this is shown for 
the version of /inorm with the 2-norm for ||/|| , which is bounded above by the expression ([T|)). 
Since Df{eo) = Dg{eo) and \\g\\ < ||/||, we can conclude /Unorm(/, eo) > A*norm(ff, eo) > 1. 

□ 



3 The exclusion Lemma 

In this article, d{ , ) denotes the Riemannian (angular) distance in S*" (which satisfies 
< d{x,y) < TT, yx,y G 5") and for x S S",r> 0, we set B{x,r) := {y G 5" : d{y,x) < r} 
and B{x,r) := {y G 5" : d{y,x) < r}. 

The following result can be used to support an exclusion test. 
Lemma 3.1. Let f G TCd and let x,y E S"" such that d{x, y) < \/2. Then, 

\\f{x)-f{y)\\oo < WfWVB d{x,y) 
In particular, if f{x) ^ 0, there is no zero of f in i3(a;, min{||/(.-E)||oo/(|l/|| \/D), \/2})- 
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Proof. An immediate consequence of the definition of the 0{n + l)-invariant inner 
product is that Tid endowed with this inner product is a reproducing kernel Hilbert space [9l 
Prop. 2.21]. This imphes that, for aU g € Tid and x € M"+\ 

g{x)^{giX),ix'^Xr^^). (3) 

Because of orthogonal invariance, we can assume that x — eg and y = eg cos 9 + ei sin 9, 
where 9 = d{x,y). Equation ^ implies that 

Mx) - My) = (MX), {x'^Xf^) - {y^Xf^) ^ {x'^Xf^ - {y'^Xf^) 

= (/, {X) , ~ (Xo cos 9 + sin 0)* ) . 

Hence, Cauchy-Schwarz-Bunyakowsky implies: 

\h{x)- h{y)\<\ I /. 1 1 1 1 - (Xo cos 9 + Xi sin 9f^ \ \ . 

Since 

Xg* - (Xocos^ + Xisin^)'^' = X^^{1 - (cose)"**) {^^{cos9Y^-^{sin9fX^^-''X^, 



we have: 



\\X^^ - [Xa cos 9 + Xi sin 9Y^f = (1 - (cos e)''')^ + ^ j^*^'^ (cos 6l)2('''-'=) (sin ^^^fc 

= (l-(cos0)''')2 + l-(cos0)2^- 

= 2(1 - (cos 6')''') 

< 2(1 -(1-^)'^') (4) 

< 2(l-(l-rf,^)) (5) 

where the inequality in line ^ is obtained from Taylor expanding cos0 around 0, and the 
inequality in line ([5]) is due to the fact that (1 — aY >1 — da iox a <1. 
We conclude that 

mx)-h{y)\ < \\m9y/d; 

and hence 

ll/(a:)-/(y)lioo < 11/11 ^?^maxd,. 

For the second assertion, we have 

ll/(y)lloo > ll/(x)|U-||/(x)-/(2/)|U 

> ||/(a;)||oo- 11/11 VDd(x,y) since ^(x, zj) < \/2 

> ||/(a;)||oo-||/||x/D||/(.T)||oo/(||/||VD) = 0. 

□ 
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4 The proximity Theorem 

4.1 Newton and Smale 

Newton iteration on the sphere S*" is defined by 



TV/ : S'" -> 5" 

a; h-^ iV/(a;) exp^ (^-i:>/(a;)|j?_5„/(x) 



where exp^ is the exponential map at x, 



exp^ h — cos{\\h\\)x - 



\h\\ 



Furthermore, the standard invariants of a-theory, introduced by Smale in '24], can be 
defined as: 



(3{f,x) = Df(x)-l^J{x) 



sup 

fe>2 



Df{x)-l^^D^f{x)\(T.S^ 



k\ 



i/(fc-i) 



a{Lx) = (3{f,x)j{f,x). 

Remark 4.1. 

(i) It is easy to see that /3(/, x) — d{x, Nf{x)). 

(ii) We will not use Newton's method in our algorithm. We are instead interested in its 

alpha theory which guarantees existence of zeros near points x with «(/, x) small 
enough. 

(iii) The Newton iteration presented above is not the iteration known as 'projective New- 
ton'. There is an alpha theory for that method, available in 

Here we use slight modifications of the quantities a, /? and 7, more adapted to our 
purposes. We set 



ll/(^)l|c 



J)3/2 

a{f,x) := (3{f,x)'j{f,x). 
The definition of 7 is motivated by the estimate of 7 2, Theorem 2 p. 267]. 

l{f,x) < l{f,x). 



which yields the lower bound 



k(/)> max 2D-3/2^(/,C). 
CI/(?)=o '^'^^ 



(6) 



We also observe that j{f,x) > since ^norm(/, 2:) > 1 and that P{f,x) < /?(/, x) 



since 



Therefore a{f,x) < a{f,x). 
4.2 Proximity and unicity from data at a point 

Definition 4.2. We say that x G S"" is an approximate zero for / if and only if the Newton 
sequence {xfej^gN, where xq :~ x and x^+i '■= Nf{xk), is defined for all k and moreover 



d{xk,Xk+i) < y-j d{xo,xi). 

The limit point ^ = limfc^oo Xk is a fixed point for Newton iteration and a zero of /. It is 
called the associated zero to x. 

In what follows we denote a := X]fc>o ~ 1.632843018 . . . and we set 

Bf{x) ■.^{yeS^\d{x,y)<aMf,x)}. 
The main technical tool in our algorithm is provided by the following result. 

Tiieorem 4.3. There exists an universal constant a* :~ 0.0384629388 . . . such that for all 
X G S", if a{f,x) < a*, then 

(i) X is an approximate zero of f. 

(ii) // C denotes its associated zero then G Bf{x). 

(iii) Furthermore, for each point z in Bf{x) the Newton sequence starting at z converges 
to C- 

4.3 Background material 



Theorem 14.31 is a consequence of the following two results, which are restatements of results 
proved in ^U\. While fTOj deals with Newton iteration on arbitrary complete real analytic 
Riemannian manifolds, here we reword them in terms of Newton iteration on the unit 
sphere 5" (Example 1 in TU]). The 7-Theorem for mappings (TUl Theorem 1.3] becomes 
the following. 

Theorem 4.4. Let f : 5" R" be analytic. Suppose that f{() = and Df{() is an 
isomorphism. Let 

i?(/,C):=min|.,|-^| . 

If d{x,C,) < R{f,C), then the Newton sequence Xk — N^{x) is defined for all k > and 
d{xk,C) ^ (5)^ ^ d(x,(). In particular, {x^} converges to C,. 
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Now let ao := 0.130716944 . . . denote the smallest positive root of the polynomial ip{u)^— 
2u, and 

So := r = 0.103621842 . 



V'(crcto) \ 1 — crao 

We state the a-Theorem for mappings 10, Theorem 1.4] for the sphere 5". 

Theorem 4.5. Let f : S" ^ M" be analytic. Let x E S" be such that P{f,x) < sqtt and 
a;) < uq. Then the Newton sequence Xk = N^{x) is defined for all k > and converges 
to a zero Co//. Moreover, 

s 9* — 1 
-^X 2 1 



d{xk,Xk+i) < y-j P{f,x) 
and 

d{xkX)<^P{f,x). 

Finally we introduce i]:{u) := l—4:U+2u^, which is positive and decreasing for < u < 1 — 
and state [10, Lemma 4.3]: 

Lemma 4.6. Let x,y E S" with d{x,y) < tt. Suppose that Df{x) is nonsingular and 

V := d{x, y)j{f, x) < 1 - 

Then 

4.4 Proof of Theorem 14.31 

Set ly^, :— 0.0628039411 ... for the only real root of the polynomial 

(3 -V7){1- u)^P{u) - Au, 

and = 0.0384629388 .... Note that a* < min{ao, sqtt}. 

Since j{f, x) > the hypothesis of Theorem [43] hold from a{f, x) < a{f, x) < a* < 

ao and /?(/, x) < P{f, x) < < < so^- 

Using Remark 14. iT i) it follows that a; is an approximate zero of /, and that the associated 
zero C satisfies: 

dixX)<aPif,x)<a^{f,x). 



This already proves Parts (i) and (ii) of Theorem 14.31 
We show (iii). Since d{x,() < crP{f,x) < asQir < tt, 

— V2 

ly = d(x, Chif, x) < aPif, x)7(/, x) < aa{f, x) < era* = i/, < 1 - — , 



and we can apply Lemma 14.61 Therefore 

W(/,.)7(/,C) < 4-/3(/,-)7(/,.)(^3^ ^ 4-- (l^..W) - 3 - V7, 
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because {l — u)ip{u) decreases for < u < 1 — ^, and i^^ is a zero of {3 — — u)^{u)~ iu. 
This shows, since 2af3{f, x) < tt, that 

2a/3(/, x) < Rif, C) = min jvi, | . 

We conclude applying Theorem 14.41 to z £ Bf{x), since 

d{z, C) < d{z, x) + d{x, C) < 2a^(/, x) < i?(/, C). 

It follows that the Newton sequence {z^jfegN starting at z converges to C- 

Remark 4.7. The hypothesis on the radius of injectivity in [10^ was recently found to be 
redundant. 

5 Infinite precision 
5.1 Grids and Graphs 

Our algorithm works on a grid on S*". We easily construct one by projecting onto S" a 
grid on the cube C" = {y \ \\y\\oa = !}• We make use of the (easy to compute) bijections 
^ : ^ and : 5" given by q>{y) = ^ and ^-\x) = pjj-. 

Given rj := 2^*^' for some k > 1, we consider the uniform grid Utj of mesh 77 on C". This is 
the set of points in C" whose coordinates are of the form i2~'' for i G {—2*^, — 2*^' + 1, . . . , 2'''}, 
with at least one coordinate equal to 1 or —1. We denote by Gi] its image by cj) in S". Note 
that, for y 1,2/2 G C", 

rfWyi)X2/2)) < |||yi-2/2||2 < |Vn + i||yi-2;2||oo. (7) 

Given 77 as above we associate to it a graph Grj as follows. We set A{f) := {x S 5" | 
a{f, x) < a*}. The vertices of the graph are the points in Qjj n A{f). Two vertices x,y £ Qjj 
are joined by an edge if and only if Bf{x) (1 Bf{y) ^ 0. 

Note that as a simple consequence of Theorem 14.31 we obtain the following lemma. 

Lemma 5.1. 

(i) For each x G A{f) there exists € Z{f) such that (^^ G B f{x). Moreover for each point 

z in Bf{x), the Newton sequence starting at z converges to (^x- 

(ii) Letx,yeA{f). Then ^ Cy ^ B f{x) DB f{y) ^ (d. □ 

We define Z{Gn) ■— [JxeG ^ fi^) "^'^ where x G has to be understood as x running 
over all the vertices of . Similarly, for a connected component U of G,j , we define 

ZiU) := U Bf{x). 

Lemma 5.2. 

(i) For each component U of Grj, there is a unique zero Cu G Z{f) such that C,ij G Z{U). 
Moreover, (u G (IxguB f{x). 
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(ii) If U and V are different components of Gn, then Cu ^ Cv ■ 

Proof. (i) Let x G U. Since x G A{f), by Lemma [5.11 (i) there exists a zero (x of / 
in Bf{x) C Z{U). This shows the existence. For the second assertion and the uniqueness, 
assume that there exist C, and ^ zeros of / in Z{U). Let x,y G U he such that ( G Bf{x), 
and ^ G B f{y). Since U is connected, there exist xq — x, xi, . . . , Xk-i, Xk '■= y in A{f) such 
that {xi,Xi+i) is an edge of G,, for i = 0, . . . , fc — 1, that is, B f{xi){^B f{xi+i) 7^ 0. If Q and 
Ci+i are the associated zeros of Xi and Xi+i in Z{f) respectively, then by Lemma ISTlT ii) we 
have Ci — Ci+i; a-nd thus C = C S B f{x) (1 Bf{y). 

(ii) Assume (^ = G Bf{x) n Bf{y) C 2'(C/) n Z{V), then x and y are joined by an edge 
and belong to the same connected component. □ 



5.2 The infinite precision algorithm 

Count _Roots_l(/) 
let rj := 

(1) let [/i, . . . , (7r be the connected components of Gn 
if 

(i) for 1 < i < j < r 

for all Xi G Ui and all Xj G Uj, d(xi,Xj) > nri\/n + 1 

and 

(ii) for alla;Ge^\A(/), ||/(x)||oo > f r?V('^ + 1)D||/|| 
then HALT and return r/2 

else 77 := r;/2 
go to (1) 

5.3 Proof of Theorem [IIl](l-3) 

Proof of Part (1) This proof requires some arguments of convexity. We can natu- 
rally define spherical convex hulls for sets of points in i/", an open half-sphere in S*". If 
xi , . . . , G i?" we define 

SCH(xi, . . . , Xq) :— Cone(xi, . . . , Xq) n S*" 

where Cone(a;i, . . . , Xq) is the smallest convex cone with vertex at the origin and containing 
the points xi, . . . ,Xq. Alternatively, we have, 

SCH(xi, ...,xq)^\ l'"'' + • • • + ^^^^ I Ai, . . . , A, > 0, ^ A. = 1 j . 

U|Aia;i H \- >^qXq\\ ^ J 

We will use the following fact. 

Lemma 5.3. Letxi,...,Xq G iJ" C If ^\i=lB{x^,r,) ^ 0, then SCH(xi, . . . , a;,) C 

\JUBix,,r,). 

Proof. Let x G SCH(a;i, . . . , Xq) and y G Hi^i B{xi,ri). We will prove that x G B{xi,ri) 
for some i. 

U X — y, this is obvious. 

li X ^ y, let H he the half-space 

H -.^{zG M"+i : {z,y~x) <0}. 
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Since = \\y\\ = 1, we have {x + y,y — x) =0, and we note that in this case, x + y 
determines the mid-hne between x and y. Moreover, since x ^ y, we have x € H since 
{x,y — x) — {x,y) — \\x\\^ < \\x\\ \\y\\ — — 0. Therefore the half-space H is the set of 
points z in such that the Euchdean distance ||z — xH < \\z — y\\. 

On the other hand, H must contain at least one point of the set {xi, ...,Xq} since if 
this were not the case, the convex set Cone(CH(a;i, . . . , x^)) would be contained in {z : 
(z, y ^ x) > 0}, contradicting x G SCH(a;i, . . . , Xq). Let, therefore, Xi £ H. It follows that 

li^; - < \\y - 

which implies 

d{x,x^) < d{y,Xi) < Ti. □ 



We can now proceed. Assume the algorithm halts, we want to show that if r equals 
the number of connected components of G,,, then #r(/) = #Z(/)/2 = r/2. We already 
know by Lemma 15.21 that each connected component U of Gn determines uniquely a zero 
Cu e Z{f). Thus it is enough to prove that Z{f) C Z{Gy^). 

Assume that there is a zero C of / in 5" such that ^ is not in Z{G,i). Let i3oo(0^^ (C)j v) •= 
{y G I \\y ~ 0"HC)l|oo <v}^ {yi, ■ ■ ■,yq}, the set of all neighbors of </>"^(C) in U,^, and 
let Xi = 4>{yi), i — 1, . . . ,q. Clearly, is in the cone spanned by {j/i, . . . , j/g} and hence 

C e SCH(a;i, . . .,Xq). 

We claim that there exists j < q such that Xj ^ A{f). Indeed, assume this is not the 
case. We consider two cases. 

(a) All the Xi belong to the same connected component U of G,, . By Lemma 15.21 there 
exists a unique zero (u € S*" of / in Z{U) and G r]iBf{xi). We may apply Lemma [Ol 
to deduce that 

SCH(a:i,...,Xg) C [jBf{x,). 

It follows that, for some i G {1, . . . , g}, C G Bf{xi) C Z{U), contradicting that ( ^ Z{Gr^). 

(b) There exist i ^ s and 1 < « < j < such that xi G IJi and Xs G Vj- Since condition 
(i) in the algorithm is satisfied, dixi^Xs) > nrjy/ri+l. But, by 

d{xi,Xs) < ^Vn+l\\ye-ys\\oo < ^Vn + 1 {\\ye - (l)-\C)\\oo + WHC) - ys\\oo) < nr^Vn+l, 
a contradiction. 

We have thus proved the claim. Let then I < j < q he such that xj ^ A{f). Since 
condition (ii) in the algorithm is satisfied ||/(a;j")||oc > + 1)D||/|| . It follows from 

the inequality d{xj,C) < fV" + ^-nd Lemma \3A\ that ||/(C)||oo > 0, a contradiction. 

Proof of Part (2) We need a few lemmas. 
Lemma 5.4. If d ^ C2 £ Z{f) then 
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Proof. For i — 1,2, using (O and Proposition 12. 21 

= rmx^)] - 1^' — <r) — 1 ^ — <r) — ■ 

Now suppose that (i(Ci,C2) < RifXi) + C2) and choose x G S*" such that < 
R{f, Ci) a-nd (i(a;.C2) < R{f, C2)- Then Theorem 14.41 imphes that Ci = C2, a contradiction. 

□ 

Lemma 5.5. Let xi,X2 £ Gn with associated zeros Ci 7^ C2- If V ^ ^37rK(/)V»+i 

d{xi,X2) > TTTj^/n + 1. 



Proof. Assume (i(a;i,a;2) < Trr/V n + 1 . Since X2 ^ Bf{xi), d{xi,X2) > cr/3(/, xi). 
Consequently, 

d(a;i,Ci) < o-/3(/,a;i) < d(a;i,a;2) < TryyVn + 1 



and, similarly, d(a;2,C2) < T^rjy/n + 1. But then, 

, 2(3 - \/7)D-3/2 

rf(Ci,C2) <d(Ci,a;i) + d(xi,a;2) + d(a;2,C2) < SnijVn + 1 < ^ -^^ 

contradicting Lemma □ 

Lemma 5.6. Let x E such that x ^ A{f). If t] < ^^+1)10'^ K(f)'^ then ||/(x)||oo > 
fr;v/(n + l)D||/||. 

Proof. Since x ^ ^(/) we have a{f,x) > a*. We divide the proof in two cases. 



Case I. min 



m |/inorm(/, x), | 



In this case 



< TT < 



(n + l)D2K(/)2 - (n + l)D2||/||2 
which implies, since 77 < | < . 



||/(x)||^ > ^ > -W(" + 1)D|| 



Case II. min 



in |^/inorm(/, 2;), ^ — Mnorm (/j 2^) 



In this case 

a* 

77 < -. < 



(n+l)D2«;(/)2 - (72+l)D2/.„,,,„(/,a;)2 
which implies a* > r/(r7 + l)D2/i„oj.;jj(/, x)2. Also, 

a,<aif,x) = i^(/,x)M„orm(/,2:)D3/2 < _i_^„^^^(/,2;)2D3/2|l/(a;)|l, 
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Putting both inequalities together we obtain 
or yet, 

||/(a;)||oo >2rKn+l)Di/2||/|| > ^^^(^ + 1)D|I/II- □ 



We can now conclude the proof of Part (2). Assume t] < („^i)d2k(/)'^ • Then the 
hypotheses of Lemmas 15.51 and | 5 . 61 hold . The first of these lemmas ensures that condition (i) 
in the algorithm is satisfied. The second, that condition (ii) is so. Therefore, the algorithm 
halts as soon as 2{n+i)o^Kify^ < V < («+i)d^k(/)2 ■ This gives a bound of 0(ln(nDK(/))) 
for the number of iterations. Since the number of grid points considered at this iteration 
('^ ^ frt+ilD^Kf fl^ ) is at most 2(n + 1) j ^ ^he bound for the total complexity 



follows. 

Proof of Parts (3) and (5) We have already seen that the number of iterations is 
bounded by C'(ln(nDK(/))). At each of these iterations, we need to perform a number of 

computations on the (at most) 2{n+ 1) '^(/) ^ gj-y points to decide whether they 

are in A{f). These can be done independently. Then, we need to compute the number of 
connected components of G,,. This can be done (see, e.g., [15]) in parallel time 0(ln(|T^|))-^ 
where \V,j\ denotes the number of vertices of G,, and therefore, in parallel time at most 
C'(n^(ln(nDK(/))^ + ln(Q!*)^)). Since this is the dominant step in the computation at a 
given iteration, it follows that the total parallel time consumed by the algorithm is at most 
C'(n2ln(nDK(/))(ln(nDK(/))2 + ln(a*)2)). This shows part (3). For part (5), just note 
that, for i = 1, . . . ,r, any vertex Xi of Ui is an approximate zero of the only zero of / in 

zm. □ 



6 Finite Precision 

6.1 Making room to allow errors 

Our finite precision algorithm will be a variation of Algorithm Count_Roots_l. But since 
finite precision computations will be affected by errors, we need to make room in the infinite 
precision algorithm to allow them. For this aim, we state the corresponding version of 
Theorem EH 

Theorem 6.1. There exist a universal constant a, = 0.028268 • • • such that, for all x G S" , 
if a{f,x) < a,, then 

(i) X is an approximate zero of f. 

(ii) // C denotes its associated zero then C G Bf{x). 

(iii) Furthermore, for each point z s.t. d{x,z) < 2aP{f,x) the Newton sequence starting at 
z converges to 
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Proof. Parts (i) and (ii) follow from Theorem 14.31 and the fact that a, < a*. Part (iii) 
is proved by taking i/, ~ 0.046158 • • • to be the only real root of the polynomial '^(u) :— 
(3 - v^)(l - u)i){u) - 6m, and a, = ^ ^ 0.028268. Then, one proves as in Theorem O 
that 3cr/9(/, x) < R{f, C) from which it follows that, for all z s.t. d{x, z) < 2a(3{f, x), 

d(z, C) < diz, x) + d{x, C) < SaPif, x) < Rif, C) 
and hence, that the Newton sequence {zk}ki£N starting at z converges to ^. □ 
The proofs of Lemmas 15.51 and 15.61 yield, mutatis mutandis, the following results. 

Lemma 6.2. Let xi,X2 G G,, with associated zeros and £,2, £,1 ^ ■^2- If V < 3^k(/)V»+i 
then d{xi,X2) > 2nrj^/n + 1. □ 

Lemma 6.3. Let x £ such thata{f,x) > If rj < 4D2(„"')^(j)2 then \\f{x)\\ao > 
TT7j^in + l)-D\\f\\. □ 

6.2 Basic facts 

We recall the basics of a floating-point arithmetic which idealizes the usual IEEE standard 
arithmetic. This system is defined by a set F C M containing (the floating-point numbers), 
a transformation r : R — > F (the rounding map), and a constant u G M (the round- off unit) 
satisfying < u < 1. The properties we require for such a system are the following: 

(i) For any a; G F, r{x) = x. In particular, r(0) = 0. 

(ii) For any a; G M, r{x) = x{l + S) with \d\ < u. 

We also define on F arithmetic operations following the classical scheme 

xoy — r(x o y) 
for any x,y £¥ and o g {+, — , x , /}, so that 

o : F X F ^ F. 

The following is an immediate consequence of property (ii) above. 

Proposition 6.4. For any x,y £ ¥ we have 

x°y — {x o y)(^l -\- S), \S\ < u. □ 

When combining many operations in floating-point arithmetic, quantities such as 
nr=i(l ~^ ^i)^' naturally appear. Our round-off analysis uses the notations and ideas in 
Chapter 3 of [17], from where we quote the following results: 

Proposition 6.5. // |i5i| < u, pi G {—1, 1}, and nu < I, then 

n 

l[{i + s,)p' = i + en, 

i=l 

where 

K\<ln = -, . □ 

1 — nu 
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Proposition 6.6. For any positive integer k such that ku < 1, let 9k,0j be any quantities 
satisfying 

1^1 ku ,^ , ju 

\Ok\<ik^z — r- \Oj\<ij- 



l-ku < ■!< - '-I 1-ju 

The following relations hold. 

1. (1 + 6lfe)(l + = 1 + Ok+j for some lOk+^l < ik+j- 
2. 

l±Ok^{ l + Ok+j ifj<k, 
l + dj \ l + 0k+2j ifj > k. 

for some \0k+j \ < Ik+j or some \0k+2j \ < lk+2j ■ 

3. Ifkuju < 1/2, then jkjj < 7min{fej}- 

4. ijk < lik- 

5. 7fc + u < 7fc+i. 

6- Ik + 1] + Iklj < Ik+j- □ 

From now on, whenever we write an expression containing 9k we mean tlrat the same 
expression is true for some 9k, with \9k\ < 7fc. 

When computing an arithmetic expression q with a round-off algorithm, errors wiU 
accumulate and we will obtain another quantity which we will denote by f 1(9)- We write 
Error (q) = \q-±\[q)\. 

An example of round-off analysis which will be useful in what follows is given in the next 
proposition, the proof of which can be found in Section 3.1 of 17J. 

Proposition 6.7. There is a round-off algorithm which, with input x,y M", computes 
the dot product of x and y. The computed value fl{{x,y)) satisfies 

= {x,y) + fi'riog^ni+ida;!, 

where \x\ = (|a;i|, . . . , |a;„|). In particular, if x = y, the algorithm computes fl(||a;|p) 
satisfying 

fli\\xf) = \\xfil + 9^,,^^„^+,). □ 

We will also have to deal with square roots and arccosinus. The following result will 
help us to do so. 



Lemma 6.8. (i) Let9 eM. such that \9\ < 1/2. Then, y/1^9 = 1-9' with \9'\ < \9\. 

(ii) Let < a < 1 and e G M such that < a + e < 1. Then, arccos(a + e) = arccos(a) -f 
V , j ^ with \v\ < \e\. 

Proof. Assume > (if < it is done similarly). By the intermediate value theorem 
we have that 1 - = 0{V^y with ^ e (1 - 9, 1). But 
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the last since f > 1/2. This proves (i). 

Part (ii) is shown similarly. Again, assume for simplicity that e > 0. Then, for some 
f e (a, a + e), 

1 V 

arccosfa + e) ~ arccos(a) — earccos'(f) — e — , = — . □ 



We assume that, besides the four basic arithmetic operations, we are allowed to compute 
square roots and arccosinus with finite precision. That is, if op denotes any of these two 
operators, we compute op such that 

5p(x) = op(a;)(H-5), \S\ < u. 

From Lemma FG.Sf i') it follows that, for all a > 0, 

\/a(l + 0k) = Va{l + Ok+i). 

Remark 6.9. Our choice of the precision u in Theorem ILir i) guarantees that fcu < 1/2 
holds whenever we encounter 9^ in what follows, and consequently, \9k\ < "fk 5: 2ku. This 
implies that in all what follows we have 7g = 0{ug) for all the expressions g we will 
encounter. 

According to the previous remark we will introduce a further notation that will consid- 
erably simplify our exposition. For all expression g, we will write 

M Oiug). 

This notation will avoid we burden ourselves with the consideration of multiplicative con- 
stants. 



6.3 The finite precision algorithm 



Our finite precision algorithm is a variation of Algorithm Count_Roots_l in Section [57 
Given x e S*" we define below fl{A'{f)) and fl{B j:(x)), which are convenient floating 
versions of the sets A'{f) = {a; £ 5" | a{f,x) < ^a,} and B'f{y) = {z G 5" | d{x,y) < 
|fT/3(/, x)} respectively. 



Given / e and x G 5", we let M e 



be a matrix representing 



and we set (Tmin(M) 

/^noim 

a{f,x) 



1 

/d7 



1 



Df{x] 



A/-i||-i. Therefore 

l/(a;)lloo 



l/ll 



nCT,ni„(Af)-^||/(x)|| 



/3(/,xKorm(/,a;)^ = \\f\\na^UMr^f{x)u'^. 
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This impUcs that 

yeB'fix) ^ d{x,y) < ^a]3{f,x) ^ <7ruin{M)d{x,y) < ^a^\\f{x)\\^, 
xeA'if) ^ a{f,x)<^ ^ ||/||n||/(a;)|UD3/2 <a.(T„,in(A/)2. 

These statements are equivalent under infinite precision, but the expressions at the right- 
hand side are more convenient to handle when working with finite precision. This motivates 
our definitions of 

fl(B'f{x)) := |y e 5" I fl(a,„i„(M)d(x,2;)) < f l(^aV^|l/(a:)|U)| 
fl(A'(/)) := {xe5"|fl(ll/ll"ll/(^)llooD3/2)<fi(a.a,„i„(Af)2)} 

We also define accordingly the graph f l(G^) whose vertices are the points in n 

f 1(A'(/)), and with two vertices x, y joined by an edge if and only if f l(i3j(x))nf (x)) ^ 
0. Its connected components are denoted by fl{U). 

Our algorithm is the following: 

Count_Roots_2(/) 
let r, - 

(1) let fl((7i), . . . ,fl(C/r) be the connected components of fl(G',,) 
if 

(i) for 1 < i < j < r 

for all Xi e fl{Ui) and all xj e fl(r/j), il{d{xi,Xj)) > f l(|7rr;^nTT) 

and 

(ii) for all xeg^\ fl(A'(/)), fl(\\f{x)\\^) > fl(^^7)V(w + l)D||/||) 
then HALT and return r/2 

else rj :— ri/2 
go to (1) 

In the rest of the section we will see that, when the precision u satisfies u < 

, this algorithm is correct and halts as soon as < 



0(D2n5/2^(J)3(log 5+„3/2D2k(/)2)) 

a, 

4D^(n+l)K(/)^ ■ 

6.4 Bounding errors for elementary computations 

The goal of this subsection is to exhibit bounds for the accumulated error in the main 
computations of Count _Roots_2. We will rely on the basic notations and results described 
in gEl 

To simplify notation, and without loss of generality, in all what follows we assume that 
11/11 = 1. We denote by 5(7id) the sphere of such systems. Also, we do not discuss in what 
follows the accumulated error in the computation of 4> : C" S". This is a minor detail 
which can be taken care of using Lemma l6.8f i). 

Proposition 6.10. Given f G S{Ti.d) and x G S"' , we can compute ||/(2^)|loo with finite 
precision u such that 

Error (||/(a;)||^) = |D + log ^1 
where S is a bound on the number of coefficients of each fi . 
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Proof. Let / — (/i, . . . , /„). For i < n write fi — J2 c-jX-' and let S be the number 
of coefficients of fi. To compute f{x) one computes each monomial cjx'^ with f\{cjx"') — 
cjx"^ {1 + 6d)- Then, one computes fi{x) to get 

fl(/,(x)) = fl(^fl(c^x'0) 

= /,; (x) + ^ CjX-'el^^ + 0iog s E \^J^' 1(1 + ^D^ ) 

where in the third line we reasoned as in the proof of Proposition 16. 71 Therefore 
Error (||/(x)||^) < | ^ cjx^^f,^) + ^?iog5 J] [c/x^Kl + C) 

< X! I'^-'l l|2;"'ll(7D +7logS +7D7logs) 

< 7D+log S 

where we used that for any x £ 5", < II E IcjI^;''!! = \\fi\\ < ll/ll — 1 and 

Proposition 16. 61 f6). The conclusion follows from Remark 16.91 □ 



Proposition 6.11. Given f e SiHd) and x G 5", let M e 



be a matrix representing 



1 



1 



Df{x) 



in some orthonormal basis ofTxS"'. Then \\M\\ < y/n. In addition, we can compute such a 
matrix M with finite precision u such that 

llError (A/)||f = |n(log5' + D + logn)]. 

Proof. Step 1: Let y = ||^_c"^^|| ■ The Householder symmetry 

Hy = In+i - 2yy* 

swaps vectors e„4_i and x, and fixes y^ . The first n columns of Hy are therefore an orthonor- 
mal basis of T,j;S"', while the last column is x. Let H e k("+i)x" denote the submatrix 
obtained from the first n columns of Hy. With that notation, we set 



M 



1 



1 



Df{x)H. 



Step 2: We claim that Pi^x ■ Ti-di ~^ ^"i fi '^Dfi{x)\T^S" is an orthogonal projection, 
in the sense that for any fixed x, the map {Pi,x)\kQr{Pi ^)-^ is an isometry. 
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Wc use an orthogonal invariancc argument. The special orthogonal group SO{n + 1) 
acts on and on isometrically as follows: to a given Q G SO{n + 1), we associate 

respectively the following isometries: 

xi-^ Qx , fi^ fio QK 

We set y = Qx and Qi = fi o . Differentiating the equality gi(Qx) = fi{x), we obtain: 

Dgi{y)Q = Dfi{x). 

When x is fixed, we can set Q conveniently so that y = e„+i. Therefore 

Dgi{en+i)Q\T^sn = Dfi{x)\T^sr^- 
Since Q{T^S") = n„+^S" we obtain 

This means that Pi,o„+i(/i ° Q*) = Pi,x{.fi)- Thus, in order to prove our claim, it is enough 
to show that /i,e„+i is an orthogonal projection. 

Since for g = Y.jgjX\ ^(e„+i) = g(^^^+(d-i)e^+^) and since Tg^+iS'" = (ei, . . . ,e„), 
we have that for any gi G Hdi , 

-Pi,e„+i(ffO = ■^^(ffi(ei + (di-l)e„+i)> ■ • ■ ,5i(e,. + (di-l)e„+i)) • 

Hence, for any g^ G ker(Pi^e„+i)"'") i-e. such that gij = for all J ^ ej + [di — l)e„+i, 
1 < j < ^! we have 

J \j) 

We conclude that Pi^x is an orthogonal projection. 

Step 3: From the previous step, for any fi G Ti-di, using the orthogonal decomposition 
fi = fi+ ft with /° G ker Pi^x and G ker P^^, we have 

ll^i,x(/i)||i = ||Pi,x(/^)||Hll/^f <ll/ilP- 

It is now immediate from Step 1 and from the definition of ||/|| = maxj that the 
Probenius norm ||M||i? of the matrix M satisfies 

n n 

\\M\\i = J2 mAfi)\\l < E 11/^11' ^ "11/11' = " 

i=l i=l 

and hence its spectral norm ||M|| satisfies ||M || < ||M||ir < ^/n. This bound is independent 

of the choice of the basis for the space T^S"^ . 

Step 4: We next present the algorithm to compute M, given / and x. This is a non- 
optimal algorithm, and can be significantly improved if more is known on the structure of 
the polynomial system /. 

We can compute each entry to^ of the matrix M as the scalar product of -^Dfi(x) 
and the jth. column Hj := {hkj)i<k<n+i of H. 
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Proceeding as in the proof of Proposition 16. 101 we can compute wti^i^) with 



On the other hand, the vector y — yi^Ig"^ j can be computed using 2n + 4 operations, and 
clearly Error (j/j) — |log(n)] for all j. Hence, for all coefficients hkj of i?. 

Error (/ife,) = Ilog(n)l. 

Applying Proposition 16. 71 we conclude 



Error {rriij ) — |D + log S + log n\ 
= [D + log^ + lognl. 



--Dh{x)) 



The second equality holds because \\Hj\\ ~ 1 since H is unitary, and because, as in the proof 
of Step 2, 





2 


\Jdi 





1 



:-D5i(e„+i) 



ll(.9i(oi+(d,-i)c„+i), • ■ ■ ,5i(dic„+i))|P < ||.9i|P < 1- 



This implies 



llError (Af)llF < |n(D + logS" + logn) 



□ 



Lemma 6.12. Let x € and M he as in Proposition \6.lR We can compute crmin(Af) = 
satisfying 

Error K,in(Af)) = [n(log5 + D + n^/^)]. 

Proof. Let E' = M - fl{M). By Proposition [6111 

ll^'ll < II^'IIf < In(log5 + D + logn)l. 

Let A4 = fl(A/). We compute crmin(A^) = ||A/^^||^^ using a backward stable algorithm 
(e.g., QR factorization). Then the computed f l((Tinin(A^)) is the exact crmin(A^ + E") for a 
matrix _B" with 

\\E"\\ < cn^u\\M\\ 
for some universal constant c (see, e.g., [HI [17]). Thus, 

f l(a^i„(M)) = fl{anun{M)) = <T,nin(X + E") = a„,i„(M + E' + E"). 

Write E = E' + E". Then, using \\AI\\ < y/E, 

\\E\\ < \\E'\\ + \\E"\\ < \\E'\\ + cnMMW < \\E'\\ + cnM\\M\\ + \\E'\\) 

= |n(logS' + D + logn)] +cn'^u{Vn+ |n(log 5" + D + log n)]]) 

= |7i(log 5" + D + log n)j + cn^u{\fn + c''un(log S" + D + n^/^)) 

= |n(log5' + D + n3/2)J 

since the hypothesis on u implies c'un(logS' + D + n'^/^) is bounded by a constant term. 
Therefore, f l(crmin(M)) = crmin(M + E) which imphes by 11, CoroUary 8.3.2]: 



Error (a,„i„(M)) < \\E\\ < [n(log 5 + D + n3/2)J. 



□ 
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Proposition 6.13. Let f G S{Tid))- Assume u < i^(^jyi,^-£, g for a small enough constant 
and let x G S*" . Then 

(i) Ifx ^ fl(A'(/)) thena{f,x) > ia.. 

(ii) //xG fl(A'(/)) thena{f,x) < a,. 
Proof. From Proposition l6.10l 

fl(n||/(x)|UD3/2) = (||/(x)||oo + ID + log51)(7iD3/2)(l + 04) 
< nD3/2||/(x)||oo + InD3/2(D + log5)l 

. Also, from Lemma [6. 121 using that crmin(-^^) < V^, 

fl(a.(T„,in(M)2) = a. (fT,,in(Af) + In(log^ + D + n3/2)])'(l + 02) 



^min 

(M)|n(log5 + D + 7i3/2)] 
> a.a„,i„(M)2 - [7i3/2(log5 + D + 713/2)]. 

Therefore, 

n||/(x)|UD3/2 + |nD3/2(D + log5)l > fl{n\\f{x)\\o.B^/')>fl{a,al,^) 

> a,c7l,^ - In3/2(log S + -D + n^/^)] 

or yet, 

n\\f(x)\\^T>'^'-a,al,^ > -(InD3/2(D + log 5)1 + [n3/2(log 5 + D + n^/^)]) 

> -In-'^D^/^ logs']]. 



l|Mnorm(/,a;), - ||/(i) 



Case I. min ■ 

In this case k(/) > j|j(;^yj| — 'And, therefore, using the hypothesis on u and the inequality 
A^(/)>1, 

[n^D^/^logSl = uO(n3D^/2log5)<if^^^^^^i^ 

K(/)n"^Diogii 

< ifO(l)n|l/(x)|UD3/2<I^L&^ 

the last by choosing K small enough. Hence, n||/(a;)||ooD'^^^ — Q^»'^min — 
which implies ^n\\f{x)\\ao'D^^'^ > a,o-i„in(A^)^, i.e., a{f,x) > 



Case II. min | (/,a;), 117(7)1^1 = Mnorm(/, 



In this case n{f) > /inorm(/, a;) = ^ ^^,f ) ■ -l^y hypothesis on w, 

In^D^/^logS] = uO{n'-D'/'logS)<K^^f^f^^^ 

K(j) n^lJ log D 
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the last by choosing K small enough. This implies n|l/(a;)|looD'^/^ — a,CTmin(A^^)^ > 
_£^£2^^(M)_ equivalently, a{f,x) > 

This shows part (i). For part (ii), one shows as above that 

n\\f{x)\\^T>^/^-a,<jl,,,<ln^B^/^logSl 
Then, one proceeds as well by considering the two cases mm ^finormif , x) , jjj^ — | = 

\\f{x)\\^ ^'^'^ |/"norm(/, a;), } ~ Mnorm (/, 2^) ■ 



□ 



Lemma 6.14. Let 2/1,2/2 G o^nd let xi — 4>{yi), i = 1,2. Then d{xi,X2) > ^.Jli+i • 

Proof. The distance d(x\, X2) is minimized at j/i = (1, . . . , 1, 1) and 2/2 = (1, . . . , 1, 1 — ?/)■ 
Let iV = n + 1. Then 

xx2 {yi,y2f 

cos{d{xi,X2)) = ■^-||2n — 

{N- rj f 



N{N -2ri + ri^) 

{N - 



iV2 - 2Nri + Nrf 



Hence 



d{xi,X2) > arccos 1 — rj'^ ^N'^ ~ arcsin ^-^ViV^~T^ > 

□ 

Lemma 6.15. Let u < „ for a small enough constant K. For xi,X2 G Gq we can 
compute d{xi,X2) such that 

Error (a(xi , a;2jj < . 

1 V 1 

Proof. Let yi = (j)~^{xi), i — 1,2, and a — cos((i(xi, X2)), i.e., 

bl|lb2|| ' 

We have, using Proposition 16.71 

fl((2/l,y2>) = (2/1,2/2) + 6'logn||2/l||||y2|| 

and f l(||2/i|j||2/2||) = ||2/i||||2/2||(l + ^'iogn)- Using now Propositions [01l6.5[ and l6.6[ it follows 
that f 1(a) — a + e with e — |logri]]. 

2 

By choosing K sufficiently small, e < i2{n+iy ■ ^l^o, from the proof of Lemma [6.141, 



a = cos{d{xi, X2)) !i A/ 1 



[n + lf 
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and hence, using that 'Jz + y < \J z + ?>y whenever < z, y < 1, we obtain 



a + e < Wl - 



rf-n 



(n + l)2 12(n+l)2 
Using Lemma H^ii) it follows that, 

arccos(a + e) = arccos(a) 



4(?1 + 1)2 



3(n+ 1)' 



Therefore, 



Error a;2)) < 



arccos(a) + |logri]] 
\/«logn 



VI- + 



^7 JJ 



□ 



Lemma 6.16. Let f e >5'(7id)- Assume that rj > 



8D^(n+l)K(/)^ 



and u < 



D2n5/2^(/)3(iog s+n3/2D2^(j)2) if smaZZ enough, and let x,y e Qr,. Then 

(i) If ye tl{B'f{x)) then d{x,y) < 2ap{f,x). 

(ii) Ify i fl(B'{x)) then d{x,y) > al3if,x). 



Proof. By Lemmas 16.121 and 16.151 (and using CTmin(-^i) < and the bound d{x,y) < 
^rjy/n + 1 which follows from ([7])), 

Error (CTmin(M)d(a;,?/)) = ?/)Error (CTmin(M)) + (Tinin(M)Error (d(a;, y))) 



y/nlogn 



n log n 
11 ^7 JJ 



= 77ln3/2(log^ + D + n3/2)J 
< In^^^ log S + n^B^K{ff} 
the last by the bounds on ry. Also, using Proposition l6.101 



Error ^-aV^\\f{x)\\^ ] < IV^{B + log S)}. 

Therefore, for part (i), 
3 

crrnin{M)d{x,y) - -<yVn\\fix)\\ryo 

< fl{<j^in{M)d{x,y)) - f 1 ( laV^\\f{x)\\^ ) + fn^/^ logS + n^B^niffj + [V^(D + log 5)1 



< In3/2 g ^ n^B^Kiffl + |V^(D + log 5)1 
= ln^^^\ogS + n^T>^K{ffj. 
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Case I. min 



In this case «(/) > jjj^ — and, therefore, by the hypothesis on u, 



K(/)n(log5 + n3/2D2K(/)2) 



2^(/) 

the last hne by taking K small enough. This implies that o'inin(A^)rf(a;, y) < 2(7^/^11 /(a;) || 
i.e., that d{x,y) < 2cr/3(/, x). 



Case II. min ^^J.norm{f,x), jjjij]^^ = Mnorm(/, a:^) 



In this case n{f) > /inorm(/, x) = ^ ) • hypothesis on u 

K 



in^ ' log S + n^T>^K{ff I = Oin^'"^ log S -\- n^Ti'^ K{ff) 



' D2n5/2K(/)3(log S + n3/2D2K(/)2) 



48D2(n + 1)3/2k(/)3 



^ Vnri ^ y/nd{x,y) ^ gmin(M)d(a;, ;/) 



by taking if small enough and Lemma 16.141 This implies that j(Tyain(M)d{x,y) < 
fcrVn||/(a;)||oo, i.e., that d{x,y) < 2(T(3{f,x). 

This shows part (i). Part (ii) is shown in a similar way. □ 



Lemma 6.17. Let u < with K small enough and xi,X2 G Gr/- 

(i) Iffl{d{xi,X2)) < fl{^7Tf]y/n+ 1) then d{xi,X2) < 2TTf]^/n + 1. 



(ii) If il{d{xi,X2)) > f l(|7r?7Vn+T) then d{xi,X2) > nriy/n + 1. 
Proof. By Lemma [6.151 and the hypothesis on u, we obtain 



Error {d{xi,X2)) = 



y/n log ri 



/ y^logri\ Ki] 



< o ^^^^ ^ < ^vV^, 

\ 1] / logn 2 



the last by taking K small enough. Also, Error {^n-qy/n + 1) < ^ir'qy/n + 1 73. The 
statement easily follows from these two bounds. □ 

Lemma 6.18. Let u < Kri-/nD ^^^^ ^ small enough, f e SiHd) and x e . 

(i) //fl(||/(a;)||oo) < fl(^^?7V(" + l)D) then < ^?7^(n + l)D. 

(ii) //fl(||/(a;)|U) > fl(^^?7V(" + l)D) t/je« ll/Wlloo > f Tyv^l^^ + IP- 
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Proof. For part (i), from Proposition lG.lOl 

||/(x)||oo<fl(|l/(x)||oo) + ID + l0g51. 

Also, 



^TTT? Vl'^ + 1)D > f l(^7r7, V(n + 1)D) - l^^{n^\)Vi\. 



Therefore, 



\![x)\\^-^'K,^^{n+\)Vi < f - fl{^7rrjV{n + 1)D) + [D + log 5 + r^^in + 1)D1 



= 0(D + log5 + 77v/(" + l)D) ^^^^ 



D + log S + ijVnD 



< (1-^)?7V(" + 1)D, 



the last by taking ii' sufficiently small. It follows that ||/(a;)||oo < t^vV + 1)D and hence, 
part (i) of the statement. 

Part (ii) is proved similarly. □ 



6.5 Proof of Theorem 11.11 (4): Correctness 

We will show that, if u < a{D^n^/^K{f)'>{iols+n'>/^-D^K{m) ' ^'^'^ ^^'^ algorithm halts with 
ry > 8D2(„°*)^(y)2 ; then the value r/2 returned by the algorithm is #r(/). This is a 
consequence of the floating following versions of Lemmas 15.11 and 15.21 

Lemma 6.19. Let f e 5(Hd), V > sn^nl'iMfy « < o{n^n^^^.imiols+n^/^nMm) - 

(i) For each x G fl{A'{f)) there exists G Z{f) such that (x G Bf{x). Moreover for each 

point z G f l(-By(a;)), the Newton sequence starting at z converges to (^x- 

(ii) Let x,ye fl(^'(/)). Then C, - Cy ^ fl(^'/x)) n f 1(5/(2;)) 7^ 0- 

Proof. (i) Applying Proposition I6.13r ii). x e fl{A'{f)) implies that a{f,x) < a,. 
Therefore, by Theorem 16.11 there exists (x G Z{f) such that (^x G Bf{x). Moreover, if 
z e f l(i?y (a;)), by Lemma [6.16f i). d{x, z) < 2a(3{f, x) and the Newton sequence starting at 
z converges to (x- 

(ii) If Cx — Cy: then B f{x) O B f (y) ^ which implies by Lemma I6.16r ii) that there exists 
zefl(B'jix))nfl(B'f-{y)). □ 

This immediately implies, using that Bf{x) C fl{B j{x)) by Lemma [6T6tii), the follow- 
ing corresponding floating version of Lemma 15.21 



Lemma 6.20. Let f e 5(7id), V > sD^nl'iMfr ""'^ " ^ ofp^nV^^ 



(/)3(logS+n3/2D2B(/)2)) 



(i) For each component il{U) of f l(Gjj), there is a unique zero Cu G ■^(/) such that 
Cu G Z{fl{U)). Moreover Cu e nxeiiiu)B f{x). 
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(ii) If fl(U) and fl(l^) are different components offl{G'), then C,ij ^ C,v ■ 



□ 



In order to show the correctness of Count_Roots_2, we only need to prove that Z{f ) C 
Z(f l(G'^)). This easily follows adapting the proof of Part (1) in Section [5?3l to this situation, 
making use of Lemma I6.20l and the facts that Condition (i), fl{d{xi,Xj)) > fl[^'nri^Jn + 1), 
implies that d{xi,Xj) > nriy/n + 1 fLemma I6.17f ii)) and Condition (ii), f l(||/(a;)||oo) > 
fl(^7r77v'(n + 1)D), implies that ||/(a:)||oo > f7?V(" + 1)° (Lemma ElUii)). 

6.6 Proof of Theorem 11.11 (4): Complexity 

We want to show that if 77 < 4D2(„"')^(j-|2 then Count_Roots_2(/) halts. Note that this 
means that 

8D2(n + 1)k(/)2 - 4D2(n + 1)k(/)2 

and hence, by 16.51 that it correctly returns #r(/). 

Because of the hypothesis on rj, the hypotheses of Lemmas 16. 2[ and 16.31 are satisfied. 
Let fl{U) fl(^) be different components of fl(G^), and therefore, by Lemma [6.201 

^ (y^ and for aU x G fl(C/), y € f 1(F), by Lemma d{x,y) > 2^^r|^/nT^ holds. This 
implies, by Lemma l6.17f i). that Condition (i) in Count_Roots_2 is satisfied. 

Consider now x ^ fl{A'{f)). By Proposition I6.13f i). a{f,x) > This implies, by 

Lemma [5751 that ||/(a;)||oo > t^vV ("- + l)Di which in turn, by Lemma[6TUi), ensures that 
Condition (ii) in Count_Roots_2 is satisfied. Hence, the algorithm halts. 

Aknowledgement. We are grateful to Andre Galligo for a helpful discussion. 
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